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Fresh or outside air needed for ventilating underground installations is introduced 
through shafts or tunnels. For a shaft or tunnel in continuous use, heat is transferred 
from the air to the rock in summer and from the rock to the air in the winter, thereby 
reducing the seasonal cooling and heating loads of the underground installation. The possible 
benefit of this tempering effect in reducing the size and operating cost of the air cooling and 
heating system of the installation cannot be realized in the design stage without a reliable 
means of estimating the influence of the tunnel. In this paper, equations and functions for 
computing the tempering effect are developed by mathematical analysis of the problem, 
assuming a sinusoidal variation of the outside air temperature, and are shown to be in 
substantial accordance with experimental results obtained in tests made in a small-scale 
model f unnel. 



1. Introduction 

When outside air for ventilating an underground 
chamber, or for other purposes, is drawn continuously 
through a long tunnel or airway in rock, heat inter- 
change between the rock and the flowing air con- 
siderably reduces, by the time the air reaches the 
far end, the fluctuations of temperature to which 
the entering outside air was subject as a result of 
diurnal and seasonal changes. The tempering thus 
effected by the rock in cooling summer air and heating 
winter air may have an important benefit in reducing 
the design size and operating cost of the air cooling 
and heating systems of an underground installation. 
Mine and tunnel air ventilation studies [1, 2] l have 
shown qualitatively the reduction of the fluctuations 
of temperature of outside air flowing into mines, 
but have ventured no analysis of this process for 
prediction of temperatures and heat flow quantities. 

This paper presents the results of a mathematical 
analysis of the heat transfer between the air and the 
rock, assuming that the temperature of the entering 
air varies periodically and sinusoidally. For the 
case of long tunnels in rock, diurnal and short-time 
climatic fluctuations are quickly attenuated and 
damped out in the first part of the tunnel. Con- 
sequently, the significant fluctuation is the seasonal 
or annual periodic variation in outdoor air tempera- 
ture, which is known to approximate a sinusoidal 
oscillation [3]. 

The mathematical solutions obtained are for the 
case of dry air. Efforts to take account of heat 
exchanges between the air and rock due to con- 
densation of moisture from the air, or its evaporation 
from the rock, introduce mathematical complexities 
that preclude analytic solutions. However, the 
results obtained for dry air are conservative in 
regard to the amounts ol* heat transferred between 
the air and rock, as the effect of latent heat transfer 
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is to increase the total heat exchanges to values 
greater than those given by the equations for dry air 
conditions. 

To determine the validity of the analytical 
treatment for design purposes, experiments were 
conducted on a small-scale model underground 
tunnel, and the data thus obtained were compared 
with the results predicted by the analytical method. 

2. Analysis 

To obtain a mathematical expression for the heat- 
transfer characteristics of the rock surrounding an 
air-intake tunnel, let us assume that the tunnel 
takes the form of a hollow cylinder. Let us further 
assume that the temperature T Q of the air entering 
the tunnel varies periodically about its mean annual 
value T t in accordance with the following relation- 
ship : 

(To—TiJ^CTo^^ — Ti) cos ut 

where the angle cot increases by 2w radians (360°) in 
1 year. 

For convenience in the development, we will use a 
subrogated temperature 6, defined as {T—T s ), 
where T s is the steady temperature of the rock 
remotely surrounding the tunnel. On this basis, 
the subrogated annual mean temperature Z - of the 
entering air is {T t — T s ); the amplitude of entering 
air temperature variation (T 0{max) —Ti) is designated 
as A. Accordingly, the equation above becomes: 



T — T s =e =di+A cos ut. 



(1) 



The hollow cylinder has a radius, a, but is other- 
wise unbounded. The surrounding solid has a 
conductivity K, specific heat c, and thermal dif- 
fusivity a. The air flows in the cylinder at a velocity 
V in the direction of the positive x-axis. It is 
assumed that in the surrounding rock, heat is 
conducted radially but not parallel to the axis of 
the cylinder, and the thermal properties are in- 
dependent of temperature. 
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Let M' be the mass of air in contact with unit 
surface area of the cylinder, c' the specific heat of 
the air, 6 X the air temperature at a distance x and 
time t, and h the coefficient of heat transfer between 
the air and the cylinder surface. The heat balance 
equation satisfied by the temperature of the air at 
distance x is then: 

where 6 a , x is the temperature at the surface of the 
cylinder at distance x. The heat transfer at a 
distance x along the cylinder surface is : 



-K(^f)=h{e x -e av 



■)• 



(3) 



Equations (1), (2), and (3) provide the boundary 
conditions for the differential equation for conduction 
of heat in the solid : 



*%.: 



dr 2 



1 dfl r ,s 

r dr 



1 aor. a 

'a dt 



=0. 



A closely approximate solution for the temperature of 
the air at distance x ) assuming steady periodic heat 
flow, is: 

T x -T s =B x =e- AA ' [A cos (<at-AB')+0 t cos AB'] (4) 

where A involves x, and A' and B' involve h. 

The derivation of eq (4) is given in section 8.1. 
The symbols are defined in section 8.3. The functions 
A' and B f are plotted on figure 1 for a useful range of 
the variables z(=a-yJo>/a.) and n(=ah/K). 

3. Model Underground Tunnel 

The model underground tunnel (see fig. 2) was 
cylindrical, 6 in. in diameter and 33.7 ft long, sur- 
rounded by 18 in. of 1:4 cement-sand mortar. 
Several thicknesses of J£-in. hair-felt insulation were 
placed on the outside of the mortar. Air was drawn 
through the tunnel by a blower and its rate of flow 
was measured by means of an orifice in the sheet 
metal duct connecting the tunnel and the blower. 
At the entrance to the tunnel was a sheet metal 
duct containing chilled water coils and electric strip 
heaters. The chilled water coils were used to cool 
the air to a temperature below the lowest tempera- 
ture to be encountered at the entrance to the tunnel. 
The electric strip heaters were governed by a pro- 
gram controller that modulated the heat input to 
the air so that the entering air temperature varied 
sinusoidally as a function of time. Thermocouples 
were located at the entrance and exit and at various 
locations in the tunnel, and in the surrounding 
cement mortar, as shown in figure 2. The thermal 
properties of the cement-sand mortar were deter- 
mined from samples obtained during pouring 
operations. 



4 

3 

B' 2 










jy 


/I0.0 
.-6.0 






A 






^^-"4.0 

^_^-3.0 












— 2.0 
1.5 


' 


d 


^>^^- 








1.0 









































^n = 20.0 














/I0.0 
• 6.0 

^4.0 


5 






> 






• 3.0 
^2.0 




/ 


'/SSjS 






^ 1.5 


A' 






/// s 






1.0 

-— 0.8 


2 


// y 


^ 








"^0.7 

— 0.6 

— 0.5 


P 


^^^ 

























0.6 
Z = a va»/a 



Figure 1. Values of A' and B' as functions of z(=aV«/ a ) 
and n( = ah/K) for use in eq (4-). 

Each test consisted of four cycles of temperature 
variation with a period of 16 hr and an amplitude 
of about 20 °F. Because the transients due to the 
initial nonperiodic condition were virtually elimi- 
nated in the first two cycles, the last two cycles were 
used for test results. 

4. Results 

Three tests, in which weight rate of air flow, W, 
was varied, were carried out on the model under- 
ground tunnel. The parameters for these tests are 
given in table 1 . 

Table 1. Test parameters 



Angular velocity, « 

Tunnel radius, a 

Tunnel length, L 

Thermal conductivity (mortar) , K. 

Thermal diffusivity (mortar), a 

Specific heat (air), c' 

Z = a-y/o}/a 



0.3927 radian/hr 

.25 ft 

33 7 ft 

0.747 Btu/hr-ft-°F 

.031 ftVhr 

.24 Btu/lb-°F 



Pertinent test data are given in table 2. The 
observed entering air temperatures were substan- 
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Figure 2. Diagrammatic sketch of model underground tunnel. 



tially sinusoidal and conformed very closely to the 
equations (within ±0.25 °F), of the form of eq (1), 
set forth in table 3. The observed values of air 
temperature at the exit of the tunnel (i=33.7 ft) 
were found to lag in phase those of the entering air 
by about 0.14 radian. All experimental quantities 
involved in eq (4) were determined directly, except 
the coefficient of heat transfer, h. Since 6 X is a 
maximum when (ait — A1V) equals zero, eq (4) serves 
for computing A', from which the values of h were 
obtained. The values of h thus computed are pre- 
sented in table 2 and are compared with values of 
the forced convection heat transfer coefficient cal- 
culated from an equation given by McAdams [4] 
which, with appropriate properties for the model 
tunnel and air substituted, is 



h = 0M6(W') - 8 . 



(5) 



Table 2. Data and results from tests on model underground 
tunnel 



Test No. 



Entering air temps., To, °F: 

Maximum 

Minimum 

Mean, TV - 

Amplitude, A _ 

Leaving air temps., Tl, °F: 

Maximum 

Minimum 

Mean 

Mean remote temp, of solid (mortar), T„ °F__ 

6i=Ti-T„, °F 

Air flow rate, W, lb/hr _ ___ 

Coefficient of heat transfer h, Btu/hr-ft^F: 

From eq (4) 

From eq (5) - 



1 


2 


109.4 


101.4 


65.2 


60.0 


87.3 


80.7 


22.1 


20.7 


95.0 


90.2 


75.0 


71.0 


85.0 


80.6 


83.2 


80.3 


4.2 


0.4 


336.4 


367.9 


1.76 


1.85 


1.69 


1.82 



100. 8 
60.0 
80.4 
20.4 



91.5 
71.2 
81.35 
82.2 
-1.8 
451.2 



2.22 
2.14 



Table :$. 


Applied entering and predicted leaving air 

temperatures 




Entering air (eq. 1) 


Leaving air, £=33.7 (eq. 4) 


Test 1 

Tesl 2 

Test 3 


7*0=87.4 1-22.2 cos «* 
To=80. 7+20.7 cos <at 
To=80.4+20.4 cos <d 


!Tx=85.1+10.0 cos (o;/-0.1417) 
Tl=80.5+ 9.7 cos (orf-OJ 125) 
Tl=81.3+10.2 cos M-0.1457) 



Substituting the derived values of h, and observed 
values of A, 6 U and the parameters from table 1 in 
eq (4), for L=33.7 ft, we obtain the time variation 
of the leaving air temperature as given by the equa- 
tions in table 3. When these equations were plotted, 
the curves followed the experimental data within 
±0.25 °F. 

The tempering effect of the tunnel on the air is 
shown by referring to test No. 3, table 3, where the 
amplitude of the air which was 20.4 °F at the en- 
trance, was reduced to 10.2 °F at the exit. 

Although the experimental data on air temperature 
satisfied the analytical solution of the problem, it 
was felt that further substantiating information was 
desirable. Equations for the time-variation of tem- 
peratures at various radii in the surrounding cement- 
sand mortar were therefore computed from eq (12) 
for test No. 3, as shown in table 4. Plots of these 
equations with associated experimental data appear 
in figure 3. 

Table 4. Temperatures in the cement-sand mortar for test No. 
3 at cross-section x = 22.7 ft, as predicted by eq (12) of section 
8.1 



At radius r=0.33 ft 
At radius r=0.50 ft 
At radius r=0.75 ft 
At radius r= 1.25 ft 



T XtT = 82.0+3.1 cos M-0.733) 
!T*.r=82.1+1.7cos M-1.169) 
T*.r=82.2+0.7 cos (wJ-1.815) 
Tx,r=82.2±0.1 
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TIME , hr 

Figure 3. Air and mortar temperature history for model underground tunnel (test No. 3, table 2). 
Experimental data for temperatures in mortar at various radii from center of tunnel are denoted by symbols. 



5. Discussion 

For the model tunnel, values of the coefficient of 
heat transfer were computed from eq (4) and the 
temperature amplitudes of the leaving air by the 
method outlined in section 4. A comparison of 
these coefficients with those determined from eq (5) 
shows agreement within 5 percent, which is within 
the probable accuracy of the thermal properties of 
cement-sand mortar as determined by tests on a few 
samples. 

For tunnels or shafts probable in underground 
installations, the diameters are likely to be greater 
than 3 ft. Values of the coefficient of heat transfer 
are not found in the literature for airways of this 
size and would be particularly difficult to predeter- 
mine for the rough and irregular surfaces usually 
resulting from the excavating process. However, 
reasonable values may be estimated from compro- 
mise considerations of values determined for air flow 
in small cylinders and air flow parallel to planes, 
and with an allowance for roughness of surfaces. 

Substitution of the parameters and the derived 
values of h into eq (4) gave the equations for the 



leaving air temperatures given in table 3. These 
equations showed differences less than 0.3 °F when 
plotted with the experimental data. Differences 
less than 0.5 °F were found between observed tem- 
peratures in the surrounding solid (mortar) , and those 
computed by eq (12), as shown in figure 3. 

Because the experimental data agree closely with 
values computed by eqs (4) and (12), the analytic 
approach developed in this paper should be appli- 
cable for determining temperatures in tunnels of 
larger size. 

6. Application to a Practical Problem 

To design heating and air conditioning equipment 
for an underground installation, it is necessary to 
predict the maximum and minimum temperatures 
of the air when it is delivered at the installation. 
From eq (4), the maximum and minimum tempera- 
tures occur when ut=AB' and cot=ir+AB / , respec- 
tively. Using eq (4), these extreme values may be 
written 



(fiL)r=*e~ AA ' (±A+0, cos AB'). 



(6) 
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Figure 1 is a plot of A' and B' as functions of n and 
z computed for a range of values expected in a 
practical problem. 

For a sample hypothetical problem, consider a 
square tunnel with nominal dimensions 7 ft by 7 
ft by 2000 ft long, continuously delivering air to an 
underground installation at a rate of 50,000 cubic 
feet per minute, or 228,000 pounds per hour. For 
noncircular cross sections, the equivalent cylindrical 
radius, a, is equal to 2 (cross sectional area) /(per- 
imeter). Assumed parameters for substitution in 
eq (6) are given in table 5. 

Table 5. Parameters for hypothetical problem, 



Thermal conductivity of surrounding rock. 
Thermal diffusivity of surrounding rock... 


iC=1.2Btu/hr-ft-°F 
a =0.035 ftVhr 
co=0.000717 radian/hr 


Coefficient of heat transfer 


h =3.43 Btu/hr-ft 2 -°F 


Equivalent radius of tunnel . 


a=3.5 ft 


Amplitude of outside air. _ _. .. 


A=45 °F 


Mean temperature of outside air. _ . _ 


Ti=52 °F 


Remote rock temperature 

Specific heat of air 


7\=55 °F 
c'=0.24 Btu/lb-°F 




W = 228,000 lb/hr 


di . . 


Ti-T s = -6 °F 


n . - - - 


ahj K= 10 




aVW«=0.5 


A' (function of??- and z) figure 1 


= 4.48 


B' (function of n and z) figure 1 . 


= 2.15 


A . 


KLjc' W =0.0439 







The extreme values predicted by eq (6) are 

{0 L ) m = + M and -39 °F. 

Because 6 is a subrogated temperature based on a 
zero value equal to the remote temperature of the 
surrounding rock, the maximum and minimum air 
temperatures at the end of the tunnel (L=2000 ft) 
are 55+34-89 °F and 55-39 = 16 °F, respectively. 
The maximum and minimum outside air tempera- 
tures (eq (1)) are 97 °F and 7 °F, respectively. 
The maximum and minimum air temperatures at the 
end of the tunnel lag those at the entrance by an 
angle equal to AB' } which in this case equals 0.094 
radian, equivalent to a time-lag of 5.5 days. 

Although the rates of heat loss and gain of the air 
in a tunnel were not discussed in the experimental 
portion of this paper, they are of considerable interest 
in the application to practical problems. For the 
hypothetical problem, the average rate of cooling 
by the rock during the cooling season (eq (14), sec- 
tion 8) is 251,000 Btu/hr, or an equivalent of about 
21 tons of refrigeration. The average rate of heating 
by the rock in the heating season is 309,000 Btu/hr, 
equivalent to the heating effected by burning about 
3.35 gallons of fuel oil per hour. The maximum 
cooling and heating rates are substantially greater 
than the average, as indicated by eq (13). 

7. Conclusions 

Mathematical equations are derived for computing 
the tempering effect of a long tunnel in reducing the 
amplitude of climatic temperature variation of out- 
side air flowing through a tunnel. Results of experi- 
ments support the application of the developed 
equations for design purn^cod. 



The functions necessary for design use of the 
equations have been evaluated for a useful range of 
the parameters and are presented graphically. 

In the design for heating and air conditioning sys- 
tems, fresh or outside air necessary for ventilation of 
an underground installation represents a portion of 
the heating and cooling load for the installation. 
Without a means for determining the temperature 
of the air at the end of its passage through tunnels 
or shafts, the designer is forced to rely upon the 
outside air temperatures for design conditions. The 
hypothetical problem of section 6 shows that the 
heating and cooling loads for design purposes are 
considerably reduced when the temperature of air at 
the end of a tunnel is predicted by the derived equa- 
tions. These equations are included in a design 
manual for military installations [5]. 

An assumption has been made which develops a 
degree of uncertainty toward the use of the equations. 
An effective coefficient of heat transfer value "h" 
suitable for a particular application may, in most 
cases, only be determined as an educated guess, be- 
cause no predictable values for large size tunnel 
airways are available from literature. 

8. Appendix 

8.1. Solution of Heat Transfer Problem 

With the boundary conditions as set forth in eqs 
(1), (2), and (3) of section 2, eq (4) may be obtained 
by assuming that the temperature in the surrounding 
solid at a cross section at a distance x from the en- 
trance is 

B rtX =Fe i(at , where F=j(r) (7) 

and also that the temperature of tin 4 flowing a ir is: 

B x =Ge i<a \ where G=f(x). (8) 

By substituting the appropriate derivatives of eq (7) 
in the differential equation lor conduction of heat 
radially in a cylindrical system having angular 
symmetry : 






r dr 



dt 



the latter becomes 



F"+-F'——=0 

r a 



This is the modified Bessel equation of zero order, 
for which the general solution given in terms of the 
modified Bessel functions I , K , is 

For F to be finite as r approaches infinity, D=0, 
and a particular solution is : 
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and from eq (7), 



F=CK,(r^) 



Differentiating, 

©,.r-^'V!«<«V!> 



(9) 



(10) 



Substituting eqs (9) and (10) in eq (3), and solving 
for C, eq (10) becomes: 



©,. 



Q h lico T r ( lioi\ 



Simplify by letting 

Z=a-ylu/a 



K (z^i)=N (z)e% (z > 



?=\- (10a) 



n=ah/K 



Substituting the above in eq (2) 



V 



dd x , dd x 



m 



dx^dt^27rM'c'a 
Substituting from eq (8), 

VQ'+iuG 



'f-r- (A' + iB') = 0. 



ff ~ (A' + iB') = 0. 



27rM'c'a 

T ... A Kx Kx 

Letting A=- — r71X/r , r = rT 7 , , 

the solution of the differential equation is : 

0=Je [-?y-A(A'+iB')~l 

Substituting in eq (9), 

d x =Je [iut-^-A(A'+iB')l 
Equation (1) states: 



Q =e t +Ae iut 



then 



P=z 



iVi(g) 

N (z) 



fl =(t> 1 (z)+~ <j> (z). 

Where values of the functions of N (z), Ni(z), ^{z), 
and 4>i(z) are given by McLachlen and Meyers [6], eq 
(10a) becomes 



\dr) r=a 



12) 



(6 X ) nPe i( ^ 

a n+Pe u » +T/ ' 



Rationalizing and further simplifying, eq (3) be- 
comes 



•K 



dd 



K6 X P—n sin /x+m cos fi 



drj r=a a 



Kd x 



dH 2 sm M 

P n 



J=A+6 i e i » t . 
Making this substitution, 

e x =e~ AA ' |~A6 z ' ( ut - < v- AB, )+e i e" i (it +AB, y\ 

Considering only the real part, 

d x =e~ AA ' |~A cos Lt-y-AB'^+di cos (^+AB / )1. 

For frequencies co under consideration, the term oox/V 
is negligible and a closely approximate solution is 

d x =e~ AA '[A cos (at-AB')+et cos AB']. (4) 

Without giving the derivation from eq (9), the 
temperature at radius, r, in the solid surrounding 
the tunnel is : 



where 



=gj (A'+iB')=h(e x -d a , x ) 



2,(P-nein„) (fi 1} 

P . n . y ° 

n+P~ 2 SU1 " 



B'-- 



2wn cos ix 

> .n _ . ' 
— ho -2 sin m 
n P 



(% 0. 



d=Ae~ AA ' 



N(R) COS ( wi +' S '— T )+f cos ( S - T ) 



V© 2+i - 2 ^ sin/i 



where 



(12) 



R=ri'u/a=rz/a 
S=*o(B)-*o(2)-A5' 

COS /U 



r=tan * 



3— sm /x 
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8.2. Air Heat Exchange 

At any time, t, the rate at which heat is lost by 
the flowing air is: 

2= WV(0 o -0s). 

The maximum rates of cooling and heating of the 
air are given by 



q m =W'c' [±A^l + e- 2AA '-2e- AA ' cos AB' 

+di{l-e- AA> cos AB')]. 
If it is assumed that cos AB'^1.0, then 

^?V(l-r^')(±A+^. (13) 

The average rate of heat flow for a given cooling 
or heating cycle is: 

([a n =W f c'(\-e- AA ') (±^+e\ (14) 

8.3. Symbols 

Unless otherwise labeled, quantities are dimen- 
sionless when consistent units are used: 



A 



A=Kx/W'c', 
, 2w(P—n sin ju) 



; plotted on figure 1 as a function 



^+--2 sin ix ' of n and 
P n 

a = radius of hollow cylinder, ft, 

= 2 (cross section) /(perimeter), for noncircular 

cylinders, 

; plotted on figure 1 as a function 



B' = 



~+--2 sin m of n and 2, 
P n 
c=specific heat of surrounding solid, Btu/lb-°F, 
c'=specific heat of air flowing in cylinder, 

Btu/lb-°F, 
e — base of Naperian (natural) logarithms, 
h= coefficient of heat transfer between air and 

surface of cylinder, Btu/hr-ft 2 -°F, 
K= thermal conductivity of surrounding solid, 
Btu/hr-ft-°F, 
M' — weight of fluid (air) per unit of surface area 
of cylinder, lb/ft 2 , 
n=ah/K, 
P=zN 1 (z)/N (z), 

#=heat flow rate from air to tunnel wall, Btu/hr, 
r= radius, ft J 



T= temperature, °F. T is temperature of air 
entering tunnel; T f is annual mean value of 
T . T s is the steady temperature of the deep 
underground rock remotely surrounding the 
tunnel, 
t = time, hr, 

V= velocity in tunnel, ft/hr, 
W' = weight rate of air flow, lb/hr, 
x= distance from entrance of tunnel, ft, 

Z = CLyjitijoL, 

a= thermal diffusivity of surrounding solid, 

ft 2 /hr, 
A= amplitude of annual variation of T about 

the mean value T u °F, 

— (-* O(max) J- i) = l/^(-*0(max) -* O(min)) j 

/x=<^i(2)+7r/4— <t>o(z), 

6= subrogated temperature, equal to (T—T s ), °F, 
0*= subrogated annual mean entering air temper- 
ature, equal to (T-T,), °F, 
co— angular velocity, radian/hr. 

Temperatures at various positions are indicated by 
the subscripts: 

0=a1 tunnel inlet (x = Q, r<^a), 
o;=a1 distance a; from tunnel inlet (r<^a), 
L=at tunnel outlet (x=L, r<^a), 
r,x=at radius r at distance x from tunnel inlet 

<r>a), 

a,x=at tunnel surface (r=a) at distance x from 

inlet . 

Tlie subscript m denotes a maximum or minimum 
value. 
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